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Abstract. Pre-thermal flow plays an important role in the flnal state evolution 
of heavy ion collisions at RHIC. We show that flow has universal features for a 
wide range of models. This significantly reduces the uncertainty in initializing 
hydrodynamic models at RHIC. 
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^ . 1. Introduction and Basic Theory 

Flow plays a pivotal role in driving soft observables at RHIC, including spectra, 
elliptic anisotropics, and two-particle correlations. For a thermalized system, flow is 
driven by the equation of state and viscosity, and thus understanding collective flow 
is synonymous with our ability to discern bulk properties of the novel matter created 
, in heavy ions at RHIC. However, significant amounts of fiow are created in the final 

stages of the reaction, t > 7 fm/c, or the pre-thermalized stage, t < 1 fm/c. For 
the latter stages, the degrees of freedom are hadronic, and the interactions become 
' largely binary, which allows one to model the latter stage rather confidently with 

hadronic cascades. However, there is significant debate regarding the microscopic 
form of matter during the pre-thermalized stage. Early flow has been shown to be 
pivotal in reproducing femtoscopic two-particle correlation data from hydrodynamic 
models [[1] [21 131 Sj. The matter might be dominated by partonic degress of freedom, 
which may be far from kinetically thermalized, to classical color fields, which carry 
no entropy. The stress-energy tensors, T^^, for these different pictures, are very 
different. The transverse pressure T^xi could be Too/3 for thermalized massless 
partons, Toa/2 for non-interacting partons, or approach Too for longitudinal color 
fields. The longitudinal pressure, Tzz, varies correspondingly, thus leading to very 
different evolutions of the energy density during the initial fm/c. In fact, for a given 
energy density Tqq, at r = 0.1 fm/c, the energy density can easily differ by a factor 
of 2 by T = 1 fm/c, depending on the state of matter during the pre-thermalized 
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stage. 

In contrast to the evolution of the energy density, the collective transverse flow, 
as defined by 

Flow^^, (1) 

varies very little for the different pictures [ [5] . Since the transverse pressures differ 
by a factor of 3 in the pictures mentioned above, this statement seems counter- 
intuitive. To understand the statement, we first note that all the pictures above 
share the property that the stress-energy tensor is traceless. During the first fm/c, 
energy densities typically surpass 10 GeV/fm'^, which is sufficiently far above Tc to 
warrant this statement, at least at the 10% level. Since T^x + Tyy + T^z = e at early 
times, and since Txx = Tyy, one can describe the initial stress-energy tensor by two 
numbers, e and k, with 

Txx — Tyy = kToo, Tzz = (1 ~ 2k)Too. (2) 

Here, k represents the transverse stiffness, with k = 1/3 for ideal hydrodynamics, 
K = 1/2 for free-streaming particles, and k = 1 for purely longitudinal electric fields, 
or purely longitudinal magnetic fields. Energy and momentum conservation give: 

a.Too - a.ToJr, (3) 

^tTqx — ^xTxx ^ ^zTxz- 

At z = 0, Tqz and Txz vanish, but for small z longitudinal collective flow {vz — z/t), 
Toz = (Too + Tzz)z/t, Txz = T^xZ I T , (4) 
which leads to the small time behavior for Tqx and Too, 

T 

drTox = -ndxToo (5) 

T 

drTooir) = --(2-2K)roo. 
r 

Combining these, one finds that the flow grows as: 

Q Tqx _ drTox Tqx drToo , , 

-too J 00 J 00 J 00 

_ Too + Tba/T ^ (2 - 2k) 

Too Too 

Defining a to describe the linear growth of the flow at small times, 

a / Tqx \ Tqx 

a = 5r — , or — = axT, (7) 
\Joo/ J 00 
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one can insert these expressions into Eq. ^ and solve for a, 



This last expression is referred to as "universal flow" . The universality comes from 
the fact that is independent of n, meaning that the initial growth of the flow does 
not depend on the transverse stiffness, and will be identical for the very different 
pictures mentioned above. However, universality is contingent on three criteria: 

1. The stress-energy tensor is traceless. This is true for non-interacting gauge 
fields, or non-interacting massless partons, and does not require thermaliza- 
tion. For the energy densities in the first fm/c, mostly above 10 GeV/fm'^, 
this is true to the order of 10%. For lower energy densities in the critical 
region, this is a poor assumption. 

2. The longitudinal motion must respect boost invariance. This should be a good 
assumption for RHIC energies and an excellent assumption at LHC energies. 
The most important part of the assumption is that the collective velocity is 
Uz — z/t, which was applied in Eq. ([¥]). 

3. The transverse stiffness k must be independent of x and y, though it can 
depend on r. Otherwise, d^T^x could not be identified as ndxToQ. 



If universality holds for sufficient time, it can be used to initialize the flow in hy- 
drodynamic models. For the Navicr-Stokes limit of viscous hydrodynamics, the 
stress-energy tensor is determined by four pieces of information: the energy density 
e, and collective velocities Ui, or equivalently Tqq and Tgi. However, more gener- 
ally, the stress-energy tensor has 10 free parameters, as the spatial components 
might differ from their Navier-Stokes values if not given sufficient time to relax to 
their correct values. This relaxation is included in Israel-Stewart approaches and is 
characterized by a relaxation time tjs [MHIMM UHl EI] ■ This relaxation time is 
determined by the viscosity and the pressure fluctuation [[T^ [T5] . 



The pressure fluctuation can be estimated by considering a quasi-particle picture 
of partons, or in principle from Lattice Gauge theory without the need to invoke 
analytic continuation. Assuming that the viscosity is within a factor of 3-4 or less 
of the KSS limit, one would expect tjs ~ 0.5 fm/c. This suggests that the Navier- 
Stokes picture should become valid at r < 1 fm/c. Thus, if universal flow can last 
all the way to r = 1 fm/c, one could confldently initialize the collective flow for the 
hydrodynamic calculation independent of the model, beginning at r = 1 fm/c. 



2. Models 




(9) 
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Fig. 1. The evolution of the electromagnetic fields and flux tubes from two oppo- 
sitely Gaussian charge densities receding at ±c for three proper times: r = 0.1 A 
(red lines), r = 0.5A (green dashed lines) and r = 2. 5 A (blue dotted lines), where 
A characterizes the transverse Gaussian size of the radially symmetric charge den- 
sities. Each quantity is plotted as a function of the position along the x axis. The 
electric field remains longitudinal and the magnetic field curve around the z axis. 
The energy density Too spreads out over time, developing into a pulse for r >> A. 
The flow Tqx/Tqq grows linearly before its magnitude saturates at unity. 
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In order to test this possibility, we compare three very different models. The first 
two models will be based on non-interacting fields, generated by charged capacitor 
plates fiying apart at ±c, with opposite charge densities ±p{x,y). From these 
currents, one can generate the electromagnetic fields, 

E,{x,y,r) - j dx'dy' p{x',y')5{(x-x'f + {y-y' f-T^), (10) 

B(x, y, r) = j dx'dy'pix', y')5{{x - x'f + {y - y'f - r^)^^^^^!^. 

For a single Gaussian charge density, p{r) = , the resulting evolution is 

displayed in Fig. [TJ Here, A is the color coherence length, or inverse saturation 
scale. For r << A, a longitudinal electric field dominates the stress-energy tensor, 
and T^^ = Tyy = Tqq, while T^z = ^Tqq. On a time scale of A, the magnetic 
field builds up, and along with it a momentum density, Tq^. — E^By. For large 
times, the field develops into a pulse of radius t and thickness A, and the effective k 
approaches 1/2. Greatly more sophisticated calculations of interacting color fields 
reveal similar evolutions of the transverse stiffness [ 14 . 

Each of the three models will have an initial energy density with the same 
transverse profile, 

^oo(x,y,r->0)^e-(^■^+'^^)/2«^ (11) 

with R — 3 fm. The first model considers a single coherent charge across the 
entire transverse area. This is clearly unphysical for a color charge, but is chosen to 
represent an extreme limit as it maintains k « 1 for the longest time. The second 
model is to consider an ensemble of point-like charges. For this case k is 1/2 for the 
entire evolution. As a third model, we consider ideal thermodynamics of a massless 
gas, and k — 1/3. 

Figure [2] shows the fiow Tq^/Toq along the x axis for three different times for 
each of the three models, and for the analytic answer of Eq. ([8|). Indeed, all 
four models are very similar, even for the larger time of 1.0 fm/c. Also shown 
are the collective velocities Ux, which are defined by the four- velocity required to 
make Tq^ vanish in the observer's frame. The velocities differ substantially between 
the models. If one were to thermalize the systems suddenly at the chosen times, 
one could find the new velocities required to maintain the energy and momentum 
densities while matching the hydrodynamic form of the stress-energy tensor. The 
velocities then suddenly change, and once again match each other, as illustrated 
in the upper panel. This emphasizes the impact of making the suddenness in the 
laboratory frame, rather than in the frame of the fluid. Otherwise, e and Ui would 
be preserved rather than Tqq and Toi- Making the change locally simultaneous 
according to a co-moving observer would insinuate an inside-out transition, with 
the inside thermalizing before the outside. The sensitivity to reference frame for 
changing k is related to the third assumption for universality, i.e., k depends only 
on T and not on x ov y. 
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Fig. 2. Results for the three models at Fig. 3. Elliptic anisotropies, defined by 
three different times, 0.3, 0.6 and 1.0 fm/cEq. (fT2|) . as a function of the proper 
as a function of the position along the x time for a source with initial Gaussian en- 
axis. The three models described in theergy density characterized by Rx = 2 and 
text are: Ideal hydrodynamics (green tri- Ry = 3 fm. Lower Panel: Sp for coherent 
angles) , and the coherent (blue squares) (blue squares) and incoherent (red circles) 
and incoherent (red circles) limits of clas- limits of non- interacting fields. The inco- 
sical fields. herent case yields = 0, exactly as one 

Lower Panel: Ideal hydrodynamics has would obtain with non-interacting parti- 
the greatest transverse radial collective cles. The solid line shows the lowest-order 
velocity even though it has the small- (in r) quadratic contribution for ideal hy- 
est transverse pressure, Txx^ of all three drodynamics. 

models. Upper Panel: Assuming that the mat- 

Middle Panel: The flow, Tqx/Tqq, is ter suddenly behaves as if it were ideal 
nearly universal for all three models. hydrodynamics at time r, represents 
Upper Panel: The collective velocity as- the anisotropy of the altered stress-energy 
suming that the matter suddenly behaves tensor, after fitting to Tqo, Tqx and T^y. 
as if it were ideal hydrodynamics at the The result is close to the quadratic form 
prescribed time. Since this is determined approximating the behavior of ideal hy- 
by fitting Tqi: /Too, it is also nearly univer- drodynamics. If one had matched to a 
sal. If the fit had been to a Navier Stokes Navier-Stokes form, rather than an ideal 
form of viscous hydrodynamics, the forms form, the three solutions would have also 
would also have converged, but to differ- converged, but to a different line, 
ent lines. 



Universal Flow in the Early Stages at RHIC 



7 



The success of matching the flow for all three very different models all the way 
to T = 1 fm/c, suggests that one can rather confidently use Eq. ^ to initialize 
the initial flow for a hydrodynamic or viscous hydrodynamic code in a way that 
is insensitive to the anisotropy of the stress-energy tensor at earlier times. For a 
viscous code, the initial collective velocity Ui would depend on the viscosity one 
chooses to apply at one fm/c, i.e., the values u[ shown in the upper panel of Fig. [2] 
would change, but they would all change to the same curve. 

Elliptic flow can also build up in the pre-thermal stage. Figure |3] shows the 
evolution of 

^ _ Jdxdy (T^^-Tyy) 
^ J dxdy {T^x + Tyy)' 

which can be identified with the elliptic flow once the system has become free- 
streaming particles. Since Sp depends on T^x and Tyy, it varies dramatically be- 
tween the three models. However, if the systems were to suddenly thermalize, and 
one were to match Too i Tq^ and Toy to the hydrodynamic form, ep changes suddenly, 
and the models jump to approximately the same value. For the case of incoherent 
color fields, which is identical to what one would obtain for free-streaming mass- 
less partons, ep is zero before thermalization. This comes from the fact that the 
contribution to the anisotropy from the flow velocity is canceled by the opposite 
anisotropy in Ty measured by a co-mover. If the system suddenly thermalizes, this 
negative asymmetry disappears and leaves a postive value for ep This demonstrates 
that the conditions for elliptic flow can be generated even while ep is zero. 

Universality has important implications for the modeling of heavy ion collisions 
at high energy. Running a viscous hydrodynamic code requires choosing an initial 
energy profile and an initial flow. The findings here show that the uncertainty of 
the flows are largely ameliorated by universality, but one should not forget that the 
initial energy density is certainly uncertain, and varies drastically for the pictures 
discussed here. For the coherent limit, the energy density stays roughly constant 
for T << R, while it fall as 1/r for the incoherent limit, and as for the 

hydrodynamic limit. Also, the compactness and anisotropics are uncertain at the 
10% level, which translate into 10% alterations of final state elliptic flow [[15] and 
correlations [31 [3] and spectra. The universality studies shown here allow one to 
state the flow at 1 fm/c with w 10% confidence, and given that the of the order of 
10% of the flow is built up during this time, one would only expect the uncertainties 
in the initial flow field to translate to « 1% uncertainty of the final flow. 
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